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Abstract 

We numerically investigate the effect of in-plane anisotropic Fermi surface (FS) on the flux- flow resistivity pf under 
rotating magnetic field on the basis of the quasiclassical Green's function method. We demonstrate that one can 
detect the phase in pairing potential of Cooper pair through the field-angular dependence of yOf even if the FS has in- 
plane anisotropy. In addition, we point out one can detect the gap-node directions irrespective of the FS anisotropy 
by measuring pf under rotating field. 
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1. Introduction 

It is important task to clarify the internal de- 
gree of freedom for the orbital part of Cooper pair 
wave function, which is the fundamental nature of 
superconductivity. We have ever proposed a new 
experimental method to detect both the phase and 
the anisotropy of pairing potential. That is, we 
have investigated the in-plane magnetic field-angle 
dependence of the quasiparticle scattering rate in- 
side a vortex core [1] and the flux- flow resistivity pf 
[2]. Through a series of our researches, we obtained 
the knowledge that the field-angle dependence of 
pi is sensitive to the phase of pairing potential 
both in the cases of an isotropic and an uniaxially 
anisotropic Fermi surface (FS). The field-angle de- 



pendence of pf has not been investigated yet for an 
in-plane anisotropic FS. 

In this paper, we investigate effects of in-plane 
FS anisotropy on the field-angle dependence of flux- 
flow resistivity /)f(aM)- Most materials have the 
anisotropy in their FS reflecting the anisotropy of 
crystal structures. Thus, it is more realistic to in- 
vestigate pf (q^m) in the case of anisotropic FS. We 
consider two model FSs with in-plane anisotropy 
and numerically calculate Pf(aM) for those FSs 
with changing the anisotropy of FS. Our numeri- 
cal results show that one can detect the gap-node 
direction by measuring pf under rotating magnetic 
field even if the FS has an anisotropy. 

2. Formulation 
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We consider a single vortex at low magnetic field 
and at low temperature. The energy dissipation due 
to the vortex flow comes from non-magnetic impu- 
rity scattering within a vortex core. Two contribu- 
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tions to the flux- flow resistivity pf are considered [2] . 
One is the quasiparticle scattering due to randomely 
distributed impurities [1,3] and the othter is the en- 
ergy scale of the quasiparticle (QP) bound states 
inside a vortex core co'o(/cf) [2,4]. Note that this en- 
ergy scale depends on the wave vector [5] within the 
framework of the quasiclassical theory of supercon- 
ductivity. 

We assume an isotropic vortex character- 
ized by the pair potential given as A(r,/cF) = 
|A(r)|d(itF)e^^(l^l).Weset|A(r)| = A^o tanh(|r|/0 
as the spatial variation of the pair potential ampli- 
tude. Aoo is the bulk amplitude, ^ is the coherence 
length, and (/>(|?"|) is the azimuthal angle in the real 
space. (i(/cF) denotes the anisotropy of the pair po- 
tential in the /c-space. kp is the Fermi wave vector. 

The expression for the flux-flow resistivity pf is 
given as [2] 

pdT) ^ -^- ^, (1) 

^Q(feF) 



where the integral on a FS with respect to 
kF is {■■■) = {l/iyo)JdSF/\vF{kF)\---. The 
area element on an anisotropic FS is dSF = 
\kF{4^k^^k)\'^ sin Okd(l)kdOk' The total density of 
state on a FS is z/q = J dSF/\vF{kF)\' The Fermi 
velocity is VfI^f) = ^k^{k)\k=kF' I^ ^^is paper, 
we use a unit system in which h = 1. Here, we as- 
sume that the system is in moderately clean regime 
and that quasiparticles with energy e = k^T pre- 
dominantly contribute to the flux-flow resistivity 
at the temperature T [6]. Within the quasiclas- 
sical theory, the momentum dependent interlevel 
spacing of the vortex bound states co'o(/cf) is ob- 
tained by using Kramer-Pesch approximation [7] 

aSC^0(/CF) = 2|d(/CF)pA^/(|/CF±||^F±(/CF)|) [2,3]. 

Here, /cf± and Vf± are the components of /cf and 
Vf projected onto the plane perpendicular to H^ re- 
spectively. The quasiparticle scattering rate inside 
a vortex core F is given by [1,3] 



r(£) = -rJl{i- sgn[d(fcF)rf(fc^)]cose) 

1 \VF±{k'^)\ M(fcF)| 

|sine||t;F±(feF)IM(feF)l 

^g-«(so,feF)g-«(So,feF) \ \ 



(2) 



and ii(s,/cF) = (2\d(kF)\/\vF^{kF)\) jl^'^ ds'/^{s') 
with A(5') = Aoo tanh(5Y^). s' is the real space 
coordinate along the QP trajectory. One can ob- 
tain further information on the expression of F in 
Ref. [1,3]. 

We consider two model FSs I, H with in- 
plane anisotropy. The model FS I is character- 
ized by the energy dispersion ei{k) = —jii — 
2t{cos{kxa) + cos{kya)} + /c^/(2m), where t and a 
are the hopping integral and the lattice constant, 
respectively, fii is the chemical potential. The dis- 
persion in the kx — ky plane is given by the the tight- 
binding (TB) model and that in the kz direction is 
free electron model. As characteristics of the TB 
model, there are Van Hove singularities in the direc- 
tion of (7r,0) and (0,7r), at which \Vke{k)\ = [8]. 
In addition to this point, the anisotropy of the FS 
grows larger gradually with increasing the chemical 
potential below the half filling. The model FS H is 
given by the anisotropic dispersion eu{k) = —/in + 
l/(2m) {kl + kl + aHkt + fc4)/2 + Sa^k^k^ + fc^} 
[9,10]. m is the mass of charge. In numerical calcu- 
lation, we set the parameter mta^ = 1. 

For an isotropic FS, the position on the FS is 
identified by the azimuthal and the polar angle 
{(j)k^Ok)' However, in anisotropic FSs, it is identified 
by 0/e, 9k and the Fermi radius |/cF(^/c,^fc)|- We can 
parametrize the Fermi wave numbers in spherical 
coordinates: kFx = |/cF(^/e, 6>/e)| cos^/e sin6>/e, /^F^y = 
\kF{(t)k,Ok)\sin(t)ksinek,kFz = |/cF(0/c,6>/e)| cos6>/e. 
Substituting these Fermi wave numbers into the 
above two dispersions eiik) and eii(/c), and using 
a bisection method, we can determine numerically 
|^F(</>fc7^/e)| such that e(k) = 0. 

For FS I, the absolute value of the Fermi velocity 
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(3) 



\vF{(t)k,Ok)\ = 2taJsiv?{kFxCi) + sin^(/cFya) + — — ^-^ttttA:^ 
and for FS H, 



FS'/ FS 



where F^ is the scattering rate in the normal state, 
0(A;f,^f) — ^v(^f) — ^v'(^f) is the scattering angle 



\vF{(t>k. Ok)\ = — |4x (1 + ^^^Ix + 3a 
m I 

+kly (1 + a'kly + ia'kl^f + fc^} yi) 

When integrating Eq. (2) numerically, we need a 
relation beween t^F±(^/c,^/e) and VF{(l)ki^k) [!]• The 
component of VF{(f)ki^k) projected onto the plane 
perpendicular to H is given by 



|^F±(0/c,^/c)| = |'^F(^fe,6>/c)| 



(a) Line-node ^|x^-/|-wave 



cos^^(/cf) 



X a/cos^ Ok + sin^ Ok siv?{(j)k — o^^^) 

\vF{(t>k,Ok)\ 



\VY^{(t)k,Ok)\ 
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(6) 
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d{kF) 
d{kY) 



- n n \ \^F{^k,Ok)\ • a • (X \(^\ 

sm6'^(/cF) = ] . , ^ x| sm6>/eSin((/)/e - Oi^)\l) 

where Oy(kY) is the angle of the QP trajectory mea- 
sured from the aM-axis. The aM - ^m plane is per- 
pendicular to H (cm 1 1 -H"). Q^M is the magnetic field 
angle measured from the (tt, 0) direction. 

In order to investigate the relation between the 
anisotropy of the pair potential and that of the 
FSs, we consider the following four model pair 
potentials: (i) Line- node 5|a,2_^2|-wave: ^(/cf) = 
I cos(2(/)/i;) sin^ Ok\-, (ii) Line-node5|a,^|-wave: d(hY) = 
|sin(2(/)fe)sin^6>/e|, (iii) d^: 
cos(2(/)/c) sin^ ^/e, and (iv) 

sin(2(/)/i;) sin^ ^fc. While s\r^i_yi\ (^i^^^ij-wave and 
dx2_y2 {dxy)-wdive pair potentials have the same 
anisotropy, only dx2_y2 (d^^^j-wave one has the 
sign change. S\x2_y2\ {dx2_y2)-wdive pair potential 
coincides with ^i^^^i {dxy)-wdive one when the pair 
potential is rotated by 7r/4 [rad] . 



3. Results 

In Figs. 1 and 2, we show the numerical results 
of the field-angle au dependece of the fiux-flow re- 
sistivity pf for the model FS I and II. The inset 
shows the schematic figure of the pair potential on 
the anisotropic FS. We fix the temperature in this 
calculation at T = 0.35Tc. Each plot corresponds to 
different chemical potential. For the energy disper- 
sion II, the chemical potential /in is fixed to /in = 9t. 

In Figs. 1(a) and 1(b), we can see that pf has the 
maximum value when H is applied parallel to the 
gap- node direction. This behavior in the s\x'^_y2\- 
wave pair is consistent with the result for both an 
isotropic and an uniaxially anisotropic FS [2]. The 
same behavior is seen also for the dx'2-y2 and the 
(ia;^-wave pair [see Fig. 2(a) and 2(b)]. However, the 
oscillation amplitude of Pf (q^m) for the dx'2-y2 and 
the dxy-^^^c pair becomes larger than that for the 
S\x2_y2\ and the S|a,^|-wave one. In addition, for the 
(i^2_^2-wave pair, a rather sharper peak appears 
when H \ \ 7r/4. These characteristics originate from 
the sign-change in the pair potential because its am- 
plitude is the same between the s\x'2-y2\ and the 
dx'2-y2 (or the s\xy\ and the dxy)-^^^c pair. 
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Fig. 1. The field angle ayi dependence of the fiux-fiow re- 
sistivity pf in the case of the (a) line-node 5|3,2_^2|-wave 
pair and the (6) line-node S|a,^|-wave one. The temperature 
is set to T = 0.35Tc. Each curve is plotted for the different 
chemical potential. For FS II, the chemical potential is set 
to /in = 9f. The vertical axis is normalized by (a) minimum 
value and (6) maximum value for each curve. 

We notice that the curve of pf (q^m) for FS I ap- 
proaches the curve for FS II with decreasing the FS 
anisotropy (/ii = — O.Olt -^ —IM) in Figs. 1 and 2. 
The anisotropy of FS II with /in = 9^ is almost the 
same as that of FS I with jii = —IM. The cusp- 
like sharp peak appearing in an isotropic FS for 
the dx2-y2-wdive pair [2] is not observed in the case 
of these in-plane anisotropic FSs. For any pairing 
states, when H \ \ gap node, the peak becomes sharp 
with increasing the anisotropy of the FS. We con- 
sider that this behavior comes from the anisotropy 
of the FS. The important point is that one can de- 
tect the gap-node direction from the field-angle de- 
pendence of the fiux-fiow resistivity even if there is 
an anisotropy of a FS. That is, a peak of Pf (q^m) ap- 
pears in the gap-node direction irrespective of a FS 
anisotropy. 

4. Summary 

We investigated the field-angle dependence of the 
fiux-fiow resistivity for two models of FSs. As a re- 
sult, we find that the maximum value of Pf (q^m) al- 
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(b) c/jcy-wave 
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Fig. 2. The field angle ayi dependence of the fiux-fiow resis- 
tivity pf in the case of the (a) <i|2,2_^2|-wave pair and the (6) 



'^\xy\ 



-wave one. The temperature is set to T = 0.35Tc. Each 



curve is plotted for the different chemical potential. For FS 
II, the chemical potential is set to /in = 9t. The vertical 
axis is normalized by (a) minimum value and (6) maximum 
value for each curve. 

ways appears when H is oriented parallel to the 
gap- node direction even if the FS has an anisotropy. 
This result is irrespective of the relation between 
the anisotopy of the pair potential and that of the 
FS [compare Fig. 1(a) with 1(b) and Fig. 2(a) with 
2(b)]. 
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